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Box No. I Basis of the report 



1 



2. 



With regard to the language, this report is based on the international application in the language in which it was 
filed, unless otherwise indicated under this item. 

□ This report is based on translations from the original language into the following language , 
which is the language of a translation furnished for the purposes of: 

□ international search (under Rules 12.3 and 23.1(b)) 

□ publication of the international application (under Rule 12.4) 

□ international preliminary examination (under Rules 55.2 and/or 55.3) 

report as "originally filed" and are not annexed to this report): 



Description, Pages 

1-15 

Claims, Numbers 

1-14 



received on 01 .07.2005 with letter of 04.05.2005 
received on 01 .07.2005 with letter of 04.05.2005 



Drawings, Sheets 

1/5-4/5 as originally filed 

5y5 received on 01 .07.2005 with letter of 04.05.2005 

□ . a sequence listing and/or any related table(s) - see Supplemental Box Relating to Sequence Listing 

3. □ The amendments have resulted in the cancellation of: 

□ the description, pages 

□ the claims, Nos. 

□ the drawings, sheets/figs ' 

□ the sequence listing (specify): 

□ any table(s) related to sequence listing (specify): 



4. 



I53 Thi^ rpnort has been established as if (some of) the amendments annexed to this report and listed^below 
had l^bS^!^^^ been considered to go beyond the disclosure as filed, as .nd.cated ,n the 
Supplemental Box (Rule 70.2(c)). 

□ the description, pages 
H the claims, Nos. 11 

□ the drawings, sheetsyfigs 

□ the sequence listing (specify): 

□ any table(s) related to sequence listing (specify): 

* If item 4 applies, some or all of these sheets may be marked "superseded." 
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Box No V R easoned statement under Article 35(2) with regard to novelty, inventive step or industrial 
applicability; citations and explanations supportin g such statement ; 

1. Statement 

Novelty (N) Yes: Claims 1-10,12-13 

No: Claims 

Inventive step (IS) Yes: Claims 1-10,12,13 

No: Claims 

Industrial applicability (IA) Yes: Claims 1-10,12,13 

No: Claims 

2. Citations and explanations (Rule 70.7): 
see separate sheet 



Box No. VII Certain defects in the international application 



The following defects in the form or contents of the international application have been noted: 
see separate sheet 



Box No. VIII Certain observations on the international application 



The following observations on the clarity of the claims, description, and drawings or on the question whether the 
claims are fully supported by the description, are made: 

see separate sheet 
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Re Item I 

Basis of the report 

1 1 The amendments filed with the letter dated 04/05/05 introduce subject-matter which 
extends beyond the content of the application as filed, contrary to Article 34(2)(b) PCT. 
The unallowable amendments concern the newly added claim 1 1 . 
Indeed it has not been found in the original description, that the method may be 
performed for " mom than one holographic images detected at the same time " and 
subsequently "superposed". • 
Although the reasoning given is the letter of reply (page 3) is considered to be correct, 
the cited passages in the description (in particular page 3, lines 30-33) do not allow the 
skilled person to im plicitly consider simultaneous (alternative^, a frame sequential 
acquisition with a single CCD and a coloured wheel may also be used) and multiple 
(alternatively, a trichrome single "image" recorded with a RGB CCD may also be 
envisaged) acquisitions. 

In addition, the wording of the claim is confusing since the method starts from 
"holographic images" to obtain "holographic image" (obtaining the "reconstructed image 
defined in claim 1 , would solve the problem). 

1 2 However from the cited passages, colour digital holography may be allowed (e.g. the 
method may be respectively performed for different wavelengths components of a 
coloured holographic image to obtain a multi-colour final reconstructed object). 



0 Reasoned statement with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

? Reference is made to the following documents: 
" D1 • KREIS T M ET AL: "METHODS OF DIGITAL HOLOGRAPHY: A COMPARISON- 
PROCEEDINGS OF THE SPIE, SPIE, BELLINGHAM, VA, US, vol. 3098, 1997, 
paqes 224-233, XP000964987 ISSN: 0277-786X 
D2- BERNHARDT M ET AL: "CODING AND BINARIZATION IN DIGITAL FRESNEL 
HOLOGRAPHY" OPTICS COMMUNICATIONS, NORTH-HOLLAND PUBLISHING 
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International application No. 



INTERNATIONAL PRELIMINARY 

SS!^ ■ PC"/0003S0 

CO. AMSTERDAM, NL, vol. 77, no. 1, 1 June 1990 (1990-06-01), pages 4-8, 
XP0001 27687 ISSN: 0030-4018 

3 1 The document D1 is regarded as being the closest prior art to the subject-matter of the 
claims, and shows (the references in parentheses applying to this document): 
a method for reconstructing holographic images (cf. abstract) wherein a digitized 
hologram is detected and digitized by an image detection device (cf. paragraph 2 and 
Fig 1c) thereby generating an array having Vr pixels (see Fig. 2) whose sizes are given 
by the sampling interval of the image detection device (cf, paragraph 2). This array is 
further processed and a discrete Fresnel transform is applied on it to reconstruct the 
object (cf. paragraph 3). 

3.2 D1 fails to teach . . 

*' firstly, that, prior to apply the Fresnel transform, the digitized hologram is 
embedded into an array having Ve values, wherein Ve is larger than Vr; and, 
* secondly, that Ve is inversely proportional to the desired pixel size to be obtained 
for the reconstructed object. 

3.3 The subject-matter of claim 1 is therefore new (Article 33(2) POT). 

3 4 The problem to be solved by the present invention may be regarded as providing a 
method of reconstructing a digital hologram allowing a better reconstructed resolution 
and a scaling effect. 

3 5 The solution to this problem proposed in claim 1 of the present application is considered 
as involving an inventive step (Article 33(3) PCT) for the following reasons: 

3 6 One the one hand, although D1 teaches in §4 to embed the digitized hologram into a 
larger zero matrix (see also Figs. 5,6) in order to perform a scaling function, the 
corresponding disclosed method is applied to digital holograms processed through the 
convolution approach, which convolution approach (detailed §4 in D1) is an mutually 
exclusive alternative to the Fresnel-approximation method (detailed §3 in D1). 
Furthermore D1 makes it clear that both methods are different, in particular with respect 
to the sizes of the "reconstruction pixels" building the "resulting image" (see abstract, 
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and equations 9,16 in D1) and also stresses that the Fresnel approach leads to s.zes 
dependency of the reconstructed objects (see abstract in D1). Therefore, it is considered 
that this teaching disclosed in D1 in the context of the convolution approach cannot 
obviously be applied/transferred to the Fresnel transform approach. 



3.7 



4. 



One the other hand, D2 teaches to generate holograms using the Fresnel-approximation 
method (cf §2 in D2) and teaches to encloses the target matrix into a larger zero array 
in order to separate the reconstructed image from the noisy first order (see Fig. 1 a, 1 b 
in D2) Since D2 generate a hologram from a target image whereas the present 
application reconstruct an object from a hologram, the scaling objective problem 
addressed in the present claim is not applicable to D2 and the skilled person would have 
no motivation to derive from D2 that the size of the larger matrix should be inversely 
proportional to the desired pixel size to be obtained for the reconstructed object 

Claims 2-10 and 12-14 are dependent on claim 1 and as such also meet the 
requirements of the PCT with respect to novelty and inventive step. 



5. The industrial applicability (Article 33(4) PCT) is clearly present for the subject matter of. 
all the claims. 



Re Item VII 

Certain defects in the international application 

6. In the original description, the following points should be checked 

* page 2, lines 6,7, M and N are inverted; 

* paqe 8, line 14, and instead of ad; 

* page 1 0, line 21 , in the equation, left hand side :delta ksi et delta eta, right hand s.de, 
indexes m,n instead of l,k. 



Re Item VIII 

Certain observations on the international application 
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7 Claim 1 does not meet the requirement of Article 6 PCT in that the subject matter for 
which protection is sought lacks clarity because the expression "equal to OS" refers to 
an unknown value "OS" which may be interpreted as Os (zeros). Correspondingly, in 
claim 2 the expression in bracket may be cancelled. 
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METHOD FOR MODIFYING SPATIAL RESOLUTION IN THE 
RECONSTRUCTION OF IMAGES IN DIGITAL HOLOGRAPHY 

5 The present invention refers to a method for modifying spatial 

resolution in the reconstruction of images in digital holography. 

More particularly, the method according to the present invention 
allows the process of reconstruction of images employed in the 
interferometric technique of digital holography to be improved, thanks to 
10 the improvement of the spatial resolution of the reconstructed complex 
field, which allows upgrading applications of Digital Holography technique. 

.The interferometric technique, allowing recording and 
•TT \Y reconstructing the complex (amplitude and phase) field reflected,. 
O transmitted and/or scattered by an object, is commonly called in scientific 

15 literature as Digital Holography, which will be hereinafter abbreviated with, 
the acronym DH (e.g. see: US Patent No. 6,262,818, to Cuche et al., 
entitled "Method for simultaneous amplitude and quantitative phase 
contrast imaging by numerical reconstruction of digital holograms", and US 
Patent No. 6,246,495, to Yamaguchi, entitled "Phase-shift digital 
20 holographic apparatus"). 

It is called digital hologram an interference pattern recorded by 
means of an integrated array of radiation detectors. 

Several methods exist allowing the numerical reconstruction of 
the complex field starting from the hologram, and in particular there are 
25 the "convolution" method and the one called the "Fresnel" method. 
r In particular, in Fresnel method, as it is known, the spatial 

^ resolution of the complex field (or also amplitude and phase) is determined 

by some parameters. Some of these parameters are determined by the 
characteristics of the integrated array of radiation detectors and, in 
30 particular, by the number of elements of which the array is composed and 
by the size of the single element. Besides, other parameters are the 
reconstruction distance, determined by distance d at which the object (or 
points of its surface and its volume), and the wavelength A of the light 
source, used for creating the hologram, which are employed in the 
35 numerical Drocess of "reconstruction. 

Commonly, in literature, the spatial resolution is quantified by 
means of the "reconstruction pixel", which is expressed as a length, and 
which will be hereinafter indicated with the acronym PR. I he dimensions 



of the bidimensional PR, A£ along the x-axis and A>/ along the y-axis. 
depend on the aforementioned parameters through the following 
mathematic formulas: 

Af. «- A,--**- < 1 > 
AT Ax Af Ay 

where Af is the number of acquired pixels (acquired by an 
image acquisition device) along the x-axis, is the number of pixels along 
the y-axis, Ax and Ay are the pixel size along the two directions of x- and 
y-axis. 

From this formula, it is clear that the complex field will have a 
PR of different value at different distances, keeping constant the other " 
parameters, and in particular the size of PR increases under the increase- 
of the reconstruction distance. In such case, the spatial resolution with 
which the complex field is reconstructed will have an inferior spatial 
resolution. On the contrary, the spatial resolution will be superior at lower 
reconstruction distance since in this case the PR size decreases. 

In other applications, as for example spectroscopic or scattering 
enquiries, it is required the recording of several holograms of the same 
object under the same conditions, but obtained with different source 
wavelengths (or with different sources at different wavelengths) (e.g. see 
the paper by M.K. Kim, "Wavelenght-scanning digital interference 
holography for optical section imaging", Optics Letters, Vol. 24, Issue 23, 
1999, page 1693). In such case, by applying the reconstruction process to 
the several holograms related to each wavelength, holograms will be 
obtained which are reconstructed with different spatial resolutions, since, 
as it appears clear from equation (1), the PR for each wavelength is 
different. In particular, the reconstruction resolution will be higher for lower 
wavelengths which give lower PR values and vice versa. 

- In the state of the art prior to the present invention, there exist 
some problems connected to the fact that the reconstruction resolution is 
rigidly determined by some parameters such as distance and wavelength. 
By way of example, some particularly problematic cases will be mentioned 
in the following. 

In some applications, digital holography is used for analysing 
variations to which the object under observation is subject because of an" 
external action (e.g. force, pressure, temperature change). I he variations 
are measured in a quantitative way by subtracting the phase maps of two 



holograms recorded with the object in two different states (for example 
before and after the external disturbance action). This technique is called 
Digital Holographic Interferometry. 

In these dynamic type observations, the distance between the 
object under observation and the detection device (e.g. a camera), at 
which the hologram is recorded, could unintentionally change, obtaining 
different holograms recorded with the object placed at different distances 
from the detection device. Hence, in order to observe the object always in 
focus it is necessary to change the value of distance to be employed in 
the reconstruction process (see the paper by Ferraro et a/, in Opt,cs 
Letters 28(14) (2003), 1257-1259) for each recorded hologram. 

From equations (1), it results that the PR value is different for ■ 
holograms reconstructed at different distances, and the spatial resolution, 
with which the object (in the complex field: amplitude and phase) is 
reconstructed in the various holograms, is consequently different. This 
avoids carrying out in a direct way a difference of phase obtained, for 
instance, with two holograms separately reconstructed at two different 
distances by actually preventing digital holographic interferometry 
technique from being applied. In fact, since reconstruction resolution is 
different in the two holograms, it is then not possible to carry out a direct 
subtraction of the phase maps (a one-to-one correspondence among the 
points of the two maps does not exist). 

In general, the change of either the wavelength or the distance 
between object and camera, may make the resolution with which it is 
possible to observe, the object worst. 

This generally prevents any direct subtraction of phase between 
the two reconstructed images for detecting and quantifying small physical- 
mechanical variations of the object. Such subtraction procedure is typically 
employed in the "holographic interferometry" technique allowing different 
states of the same object to be compared. 

Similarly, in case of applications of colour DH with use of 
different wavelengths, images reconstructed with different wavelengths do 
not overlap since the PR of each reconstruction is different (e.g. see the 
paper by l.Yamaguchi, "Phaseshifting color digital holography', Opncs 
Letters Vol. 27, Issue 13, July 2002, page 1108, and the paper by J. Kato 
et al„ '"Multicolor digital holography with an achromatic phase shifter", 
Optics Letters, Vol. 27, Issue 16, 2002, page 1403. 



M. Bernhardt et al. in the article "Coding and binarisation in 
digital Fresnel holography", Optics Communications, North-Holland 
Publishing Co. Amsterdam, voL77 no.1 , Uune 1990. p.4-8. deals with the 
use of Fresnel Transform in the reconstruction of images in DH (improving 

5 diffraction efficiency). 

T. M. Kreis et al. in the article "Methods of Digital Holography: a 
comparison", Proceedings of the SPIE, SPIE, Bellingham, VA, US, 
vol.3098, 1997, pages 224-233, ISSN: 0277-786X, teaches to embed a 
digitised hologram into an array having a larger number of pixels in the 
10 specific context of the convolution approach, in order to vary the size of 
the reconstructed images. 

r , r . | n the latter article, moreover, it is concluded that the Fresnel- 

■ approach cannot be applied if reconstruction in different depths are to be' --■ 
compared, and inversely, the convolution approach is not much suited if 
15 the whole possible field of view for opaque or transparent objects has to 

be reconstructed. 

A general method of varying the resolution of the reconstructed 

image is not provided in the prior art. 

The Applicant does not know effective solutions to the above 

20 problems. 

In fact, as it results from literature, the methods presently 
employed for obtaining a better resolution in observation of objects make 
use of complex experimental apparatuses requiring particularly delicate 
calibration procedures (e.g. see: Indebetouw et al., Appl. Phys. Lett. 75, 

25 (1999)2017-2019). 
Q It is an object of the present invention to provide a method of 

reconstruction of the holographic image starting from a digitized hologram 
solving the above drawbacks and enabling further uses. 

It is also an object of the present invention to provide 
30 apparatuses and tools necessary for the execution of the method that is 

object of the invention. 

It is further object of the present invention an apparatus for 
acquiring and reconstructing holographic images making use of the 
method that is object of the invention. 

It is specific subject matter of this invention a method for the 
reconstruction of holographic images, the holographic image being 
detected by an image detection device, the holographic image being 



35 
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transformed in a digitized hologram, the digitized hologram being 
comprised of a number V r of elementary pixels, the size of which being 
. equal to the holographic image sampling intervals, and of the V r values 
respectively associated to the elementary pixels, the method comprising a 
5 first step of processing the digitized hologram array, and a second step of 
hologram reconstruction in the observation plane starting from the 
digitized hologram processed in the first step, the method being 
characterised in that the second step is carried out through discrete 
Fresnel Transform applied on an array of V e values corresponding to pixels 
10 having size equal to that of said elementary pixels, wherein said array of 
V e values) includes said array of V r values and an integer number p = V e - 
• - -■-■< -.- ■ - v r > 0 of constant values equal to OS, said number V e of values being 
>^ r ; - inversely proportional to the desired pixel size to be obtained for the 

w reconstructed image. 

15 Preferably according to the invention, said p constant values 

are null values (OS= 0). 

..... Still preferably according to the invention, said p values are 
arranged externally to said array of Rvalues. 

Always according to the invention, said p values may be 
20 arranged in a symmetrical way or in a non-symmetrical way. 

Preferably according to the invention, the digitized hologram is 
a rectangular array of V r = N r -M r values, each value corresponding to a 
rectangular pixel of sizes Ax, Ay. 

Still preferably according to the invention, the hologram 
25 reconstructed in the second step is represented by a rectangular array of 
f Ve = Ne . Me values, each value corresponding to a rectangular pixel of sizes 

a\ = (Wax) and Arj = (MM e Ay), X being the wavelength of the wave 
beam striking the object of which the hologram is recorded, and d the 
distance between the detection device and the object of which the 
30 hologram is detected, A£ and Ar, being the reconstructed holographic 
image sampling intervals. , 
According to the invention, formulas N e = {Ml Ax 2 ), M e = (XdlAy-), 

Af = Ax, Ar} = Ay may be valid. 

Advantageously according to the invention, after the second 
35 step, if each holographic image sampling interval is not equal or less than 
a certain threshold, the number of values p added to the digitized 
hologram" array is increased and the second step is carried out again. 
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Preferably according to the invention, said threshold is a 
function of the signal-to-noise ratio of the holographic image. 

According to the invention, the method can be performed for 
more than one holographic images detected at the same time for different 
5 wavelength X, said more than one images being subsequently superposed 
in order to obtain a multi-colour final holographic image. 

It is further specific subject matter of the present invention a 
computer program characterised in that it comprises code means apt to 
execute, when running on a computer, the method subject of the 
1 0 invention. 

It is still specific subject matter of the invention a memory 
,: 7 -,,-,..,.-..,r— • medium, readable by a computer, storing a program, characterised in that" 
■ ■■■ > -*v - the program is the computer program subject of the invention. 
. ' It is further specific subject matter of the invention an apparatus 

15 for detection of holographic images, comprising a digitized hologram 
processing unit, characterised in that the processing unit processes the 
detected data by using the method subject of the invention. 

The invention will be now described, by way of illustration and 
not by way of limitation, by particularly referring to the drawings of the 

20 enclosed Figures, in which: 

- figure 1 shows a block and flow hybrid diagram of the 
traditional holographic reconstruction method; 

- figure 2 shows a block and flow hybrid diagram describing the 
holographic reconstruction method according to the invention; 

25 - figure 3a shows the effect of the reconstruction in amplitude of 

Q a Talbot effect Ronchi grating made through the traditional method ; 

- figure 3b shows the effect of the reconstruction in amplitude of 
a Talbot effect Ronchi grating made through the method according to the 
invention; 

, 30 _ figure 4a shows a particular information related to the 

reconstruction of figure 3a, in relation to a certain reconstruction distance, 
as a function of the number of pixels; 

- figure 4b shows a particular information related to the 
reconstruction of figure 3b, in relation to a certain reconstruction distance, 

35 as a function of the number of pixels; .and 

- figure 5 shows a preferred arrangement of the null pixels used 
in the method according to the invention. 
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As mentioned before, digital holography consists of recording a 
distribution of interference, which is created between an object beam (that 
has interacted with the object under observation) and a reference beam, 
through an ad hoc system for acquiring images. 
5 " Such interference distribution is processed through process.ng 

methods apt to reconstruct an image of the object under observation. 

In particular, the recorded hologram is multiplied by a digital 
replica of the reference beam and the diffraction integral of this product is 
calculated. Such hologram allows a reconstruction of the object under 
10 observation to be obtained. 

The reflection holographic recording apparatus may be for 
• instance of the Mach-Zehnder type. Once analogue data are acquired," 

('""Y:^ they are processed by a processing unit. 

Making reference to figure 1, such processing unit processes 
15 data according the traditional method. The unit of preparation of hologram 
acquisition conditions or "set-up" 2 collects radiation 4 coming from source 
- . 1 and urinate with radiation 5 the object 3 under observation. Also, .n 
such set-up 2 is present a device for creating, from the beam 6 that .s 
• reflected, transmitted or scattered by the object 3, an object beam O and 
20 a device for creating a reference beam R. The object beam O and the 
reference beam R are combined in the set-up 2 so as to create an 
interference distribution 7 in a plane. Such interference creates the 
hologram 8 of the object 3 under observation, and it may be described ,n 
terms of bidimensional distribution of intensity: 

25 

O H(x,y)=\R\ 2 +\0\ 2 + R-0 + RO 

where R* and O* represent the conjugate complex of the 
reference beam and of the object beam, respectively. 

It is now necessary to specify that, as it will be shown later, the 
30 method according to the present invention is not restricted to the optical 
field and it may be applied for the numerical reconstruction of holograms 
recorded with any type of electromagnetic (for instance X rays) and non- 
electromagnetic radiation (for instance electron beams and/or acoustic 
waves), in particular, source 1 could be also made of a combination of two 
35 or more wavelengths. For this reason, type, wavelength and conerence of 
source 1 could be any. 
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The hologram 8 is acquired, digitized and stored through an 
acquisition system 9. To this end. any type of existing or future image 
acquisition system may be used. 

The acquisition system 9 internally has a device for digitizing 
5 and computer storing the acquired image 8. The digitized image is called 
"digital hologram" 10 and it is described by an array H{nAx,mAy) of NM 
numbers, obtained by the bidimensional spatial sampling of the hologram 
H(xy) 8. 

Such bidimensional spatial sampling may be described by the 
10 following formula: 



H(nAx, mAy) = H(x, y)rect{^ , ^)Z2>(* " V ~ m ^ 



where S(x,y) is a bidimensional Dirac delta function, n and m are 
integer Ax and Ay are the sampling spacings along the x-axis and the y- 
. 15 axis respectively, (NAx)*{MAy) is the area of the image of the acqu.red 
• hologram, nctfay) is a function the value of which is 1, if the point of 
coordinates (*o0 is within the part of the acquired hologram, and 0 
otherwise. 

For a perfect reconstruction of the object image,, it is necessary 
20 that the digitization process satisfies the sampling theorem. In particular, it 
has to be satisfied the condition that the spacing between the fringes 
present in the interference distribution 7 is larger than at least two pixels of 
the acquisition system 9. Hence, the sampling theorem establishes the 
minimum resolution that is obtainable with a certain experimental set-up 2. 
Q 25 One of the great advantages offered by the digital holography is 

the fact that it is possible to directly act on the digitized hologram 10 of the 
object 3 for carrying out operations on the acquired information. 

This means that different processings of the images 1 1 may be 
made on the digitized hologram 10. Through such processings, it is for 
30 example possible to eliminate zero order diffraction present in hologram 
reconstruction, or to eliminate any "phase aberration" introduced by the 

used optical system. 

The term "phase aberration" means a deformation of the wave 
front travelling through the hologram creation ad recording system. The 
35 phase aberration correction compensates such deformations and allows 
obtaining a correct reconstruction of the observed object. 
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The process of numerical reconstruction 13 of the object under 
observation is based on two steps. In the first one. the "processed" 
digitized hologram H(n,m) 12 has to be multiplied by a digitized replica of 
the reference beam R, obtaining the following formula: 

5 

F(nAx, mAy) = H (»Ax, mAy) ■ R{nAx, mAy) = 
= R\R\ 2 +R\0f +RR'0 + RRO' 

where the first two terms correspond to the zero order 
diffraction, and from the third and/or fourth term it is possible to obtain the 
image of the observed object. 
1Q The second step of the propagation process consists of the 

—'^""propagation of the field distribution F(n,m) from the plane wherein the 
Q ' camera is placed to the observation plane. This process leads to the 

■ reconstructed image 14. 

It is then possible to numerically act on the recorded and stored 
digitized hologram through an electronic device for image acquisition 
(hereinafter generically called as camera) made of a discrete set of 
sensitive elements arranged in the shape of array of 7^ rows and M 
columns, in order to obtain a higher spatial resolution with respect to the 

techniques presently in use. 
20 In order to overcome the aforementioned drawbacks of the 

traditional method, the method according to the present invention is based 
on the extension of the array of the object hologram by introducing a 
number of additional fictitious points, the intensity of which is set to zero. 

The object is then reconstructed with the technique of the 
O 25 hologram numerical propagation from the camera plane along the distance 

separating the object from the same plane of the camera. 

The hologram propagation occurs by using the bidimensional 
Fresnel transform. The advantage of such integral is that its computation 
is simple and may be very fast performed by using a discrete formulation 
expressed in terms of Fourier transform. In fact, it is well known (see 
Goodman, "Introduction to Fourier Optics", MacGraw-Hill Companies Inc., 
2nd ed 1996) that the phenomenon of light propagation from a starting 
plane to a parallel plane placed at a distance d may be interpreted as a 
space-invariant linear system characterised by a transfer function having a 
35 finite band amplitude. Such transfer function has unitary module and 
phase depending on the spatial frequencies corresponding to the two 
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orthogonal directions within the plane placed at a distance i from the 
starting plane. 

In case of propagation of an optic field through the Fresnel 
numerical integral, the transfer function phase quadratically depends on" 
5 the spatial frequencies. Consequently, dispersive effects are introduced in 
the hologram numerical reconstruction process which increase with the 
increase of the reconstruction distance and which generally contribute to 
make the reconstructed hologram spatial resolution worst. 

As it will be clarified in the following, the extension of the 
10 dimension of the hologram array, by adding null elements, allows acting 
on the size of the minimum element composing the object reconstructed 
— ■ j m age (the "reconstruction pixel"), rather improving its resolution.' ~ 

Q Making reference to figure 2, as in the traditional case, for 

carrying out the method according to the invention it is first of all 
15 necessary to have a holographic system for creating a hologram of the 

object under observation. 

Such hologram is digitized and computer stored through a 
camera. The digitized hologram is a rectangular array obtained by 
sampling the hologram by means of the camera with a step Ax along the x- 
20 axis and a step Ay along the y-axis (Ax and Ay coincide with the camera 
pixel size) for a number of points equal to N-M (N is the number of camera 
pixels along the x-axis and M is the number of camera pixels along the y- 
axis). 

The array size related to the digitized hologram is then enlarged 
25 by adding a suitable number of points so as to obtain the desired 
resolution in the hologram reconstruction process. 

By using such extended array, it is then possible to exploit the 
technique of the bidimensional Fresnel transform for reconstructing the 
image of the object under observation, so gaining in definition. 
30 ^ In the configuration preferred by the inventors, set-up 2 ,is 

designed so as to produce a "Fresnel hologram", term indicating a 
hologram that may be reconstructed through Fresnel scalar diffraction 
approximation. 

The advantages of such approximation derives from the fact 
35 that its computation is very simple and may be performed in a very fast 
manner. In case of Fresnel approximation, numerical reconstruction of the 
hologram 12 will "be carried out according to the invention through a 
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discrete formulation of the Fresnel integral expressed in terms of discrete 
Fourier transform, that is: 

H lAx,kAy) = Ae^^^ 

5 where X is the wavelength of source 1, A is a complex constant, n, m. I, k 
are integer (-M2<n,/<A/72 and -M/2<m,k<M/2), DFT is the discrete Fourier 
transform, which may be fast computed by making use of multiple FFT 
(Fast Fourier Transform) algorithms reported in literature, Ax- and Ay are 
the sampling spacings of the hologram 12 (hence in the camera plane), d 
10 is the distance between the camera plane and the observation plane, and, 
...«i,.ov r -:. r .^— finally,- and Arj represent the sampling spatial intervals in, the- 
• » >»»*.• • . observation plane which are defined by: 

0 . A# = -^- A*-^- (D 

7/ Ac M Ay 

Hence, the reconstructed object will have size (N-AQx(M-&r\). 
15 The intervals described by equations (1) substantially define the resolution 

' of the reconstructed object 14. 

As it may be noted from the preceding formula, the resolution 
also depends on, besides the number of pixels and the resolution of the 
acquisition system 9, the wavelength X of the source 1 and the 
20 reconstruction distance d. 

In the reconstruction processes, it is generally A£> Ax and A 77 > 
Ay, i.e. the reconstructed object image is characterised by an inferior 
resolution with respect to the one with which the hologram 8 has been 
digitized and recorded. 

^ 25 As shown in the following, the method according to the 

invention allows solving the aforesaid problem and, keeping constant 
wavelength and reconstruction distance d, improving the resolution of the 
image of the reconstructed object 14. 

Always making reference to figure 2, on the one hand, the 
30 proposed method allows the system of creation and recording of the 
hologram 8 not to be modified, and, on the other hand, it is compatible for 
applications wherein the hologram recording has to be carried out in real 
time and in a continuous way. 
- In the version presently preferred by the inventors, the method 

35 acts on the processed hologram 12, that is .on the hologram which has 
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been already subject to processing for a correct reconstruction of the " 
object 3. 

The size of the array describing the digitized hologram 12 is 
expanded during step 15 by adding a certain number of points as 

5 determined in step 16. 

The number of points to be added 16 is determined by the 
resolution 17 that is desired to obtain in the reconstruction process 13. 

The value of such resolution 17 may be given either by 
conditions established by requirements 18 external to the reconstruction 
10 process 13 (for instance for observing with higher accuracy the object 
image) or for compensating the loss of resolution 19 due to the 
" — reconstruction process in applications requiring, in particular, the variation ^ 
O ' 20 of A and/or d. 

^ ■": The number of points is not the only feature allowing a 

15 resolution improvement to be subsequently obtained. 

In fact, it is needed to place the introduced fictitious points (i.e. 
the difference between the points computed as above and the acquired 
points) in a suitable way with respect to the detected points. 

It is needed to be sure that the introduced zeros do not result in 
20 a transformed image (according to the function w (lAx,kAy) described 
above) presenting false frequencies. 

For example, placing zeros among not null values of a 
sinusoidal plot would clearly introduce frequencies far apart from the one 
of the sine. 

25 Although single particular arrangements could be suitably 

O ^ adopted in specific cases, the preferred arrangement according to the 

invention is the one having the fictitious points as contour of the detected 
image, that is without interspersing them among the effective points. 

Making reference to the example of figure 5, the fictitious points 
30 50 are arranged symmetrically with respect to the contour of pixels 51 of 

the detected image. 

This contour arrangement is proper to images in any number of 

dimensions. 

The classical reconstruction process 13 is then applied to the 

35 expanded array of the hologram. 

The expansion 15 of the array size of the hologram 12 allows 
obtaining, thanks to the DFT properties, a reconstructed image with a 



lower spacing (that is with a better resolution) with respect to the 
reconstruction obtained without expanding the array. 

In other words, the traditional reconstruction process based on 
the Fresnel transform implies a degradation of the resolution with which 
the object is reconstructed; instead, the addition according to the invention 
of new elements in the hologram array allows such resolution loss to be 
corrected at most obtaining a resolution equal to the physical one 
established by the sampling theorem. 

In particular, if it is desired to obtain an image of the 
reconstructed object 14 with the same resolution of the digitized hologram 
10, keeping constant wavelength X and reconstruction distance d, it is 
necessary to expand the array of the hologram 12 from size N-M to size*' 
(dXlkx 2 \{dXIky 2 ), as it is obtained by inverting formulas (1 ). 

It would be theoretically possible to indefinitely improve the 
resolution of the holographic image by still adding new fictitious points. 
Actually, since the Fresnel transform re-distributes the intensity over all the 
points, beyond a certain number of fictitious points the intensity of many 
ones would be lowered below the background noise of the signal or the 
statistical noise of the same. However, it is matter of simple computation to 
determine the maximum number of fictitious points usable in any specific 
situation. 

In figure 3 an example of application of the aforesaid method is 
reported. Such example is related to a coherent and monochromatic light 
source 1 with emission wavelength 1=532 nm, the observed object 3 is a 
Ronchi grating with step ,4=6.25 lines/mm, the acquisition and storing 
system is a CCD with 7V=512 and M=512 pixels and with square pixel size 
equal to Ax=Aj/=6.7um. A Ronchi grating illuminated with monochromatic 
light generates the known Talbot effect, i.e. by observing the light 
scattered by the grating at increasingly long distances the grating rows 
appear increasingly defocused, save at particular distances (multiple of 
the so-called Talbot distance, i.e. A 2 /A), where the rows appear well 
defined and focused again. Hence, by using such effect, we are allowed 
not to handle focusing problems. 

In particular, the holographic reconstruction of a single line of 
the aforesaid grating, the hologram of which has been recorded at an 
effective distance of 170mm, is reported in figure 3a for several values of 
the distance d between hologram plane and observation plane. 
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The marks placed on the figure show the distances at which, 
due to Talbot effect, it is necessary to observe the grating rows well 

defined and focused. 

The typical "trumpet" shape of figure 3a, obtained by varying 
5 distance d, is an indication of the reduction of the reconstruction pixel, 
according to equations (1 ). 

The reduction of the reconstruction pixel, and hence of the 
reconstructed image resolution, under increasing distance d, prevents the 
grating rows to be sharply observed; hence, there is an information loss 
10 with the increase of the reconstruction distance. 

The method according to the invention allows such information 
tQ kg recovered. 

> : ---v- Tne holographic reconstruction of the expanded hologram of 
the Ronchi grating is reported in figure 3b for the same values of distance 

15 d used in figure 3a. 

In particular, hologram size has been increased by 512 points 
along both x-axis and y-axis, hence obtaining a hologram of 1024-1024 
pixels. Obviously, the obtained shape is still a "trumpet" one, but it is 
possible to observe that by increasing the reconstruction distance it is still 
20 possible to determine distances at which the grating rows well defined and 
focused again. For better pointing out the advantage given by the present 
reconstruction method, a line is reported in figure 4a and figure 4b at a 
certain reconstruction distance (J=434 mm), related to figure 3a and figure 
3b, respectively. The difference between such distance and the recording 
25 distance (170 mm) is multiple of the Talbot distance, and the grating rows 
O should then appear well defined. 

But observing figure 4a, it is noted that the loss of resolution 
does not allow the grating rows to be sharply distinguished. 

The application of the method subject of the present invention 
30 allows overcoming such degradation. In fact, observing figure 4b, wherein 
the reconstruction of the grating expanded hologram has been carried out, 
it is noted that the grating rows are well visible and sharp. 

The reconstruction method according- to the invention 
represents a significant improvement with respect to the other methods 
35 present in literature. In fact, the method according to the invention acts on 
' the object digitized hologram, and it may adapt the resolution of the object 
reconstructed image to the various requirements of multiple applications. 



The invention is particularly, but not exclusively, intended for 
reconstruction processes in digital holography where there is the need to 
improve the resolution with which the complex field (amplitude and phase), 
transmitted or reflected or scattered by the object, is reconstructed, or to 
5 keep constant such resolution when the variation of other parameters, of 
which the same reconstruction is function, would tend to make the 
resolving power of the holographic technique worst. In numerous, above 
all metrological, applications, there exists the need to improve the 
resolution with which an object is observed, modifying as less as possible 
10 the observation apparatus and preventing the acquisition times from 
increasing. This last requirement is particularly felt in all those applications 
j^Y"-"""' where a' real time observation of the object is required. 
Y~) ■ ' ■ ' The preferred embodiments have been above described and 

some modifications of this invention have been suggested, but it should be 
15 understood that those skilled in the art can make variations and changes, 
without so departing from the related scope of protection, as defined by 
. the following claims 



o 



15 



10/566036 

(AP20 Rgc'S ^I?T0 2 5 JAN 2006 

CLAIMS 

1. Method for the reconstruction of holographic images, the 
holographic image being detected by an image detection device (9), the 
holographic image being transformed in a digitized hologram (10), the 

5 digitized hologram (10) being comprised of a number V r of elementary 
pixels, the size of which being equal to the holographic image sampling 
intervals, and of the V r values (51) respectively associated to the 
elementary pixels, the method comprising a first step (1 1 ,12) of processing 
the digitized hologram array, and a second step (13,15,16,17,18) of 

10 hologram reconstruction in the observation plane starting from the 
digitized hologram processed in the first step, the method being 

v - characterised in that the second step is carried out through discrete 
r ~y Fresnel Transform applied on an array of V e values corresponding to pixels 

having size equal to that of said elementary pixels, wherein said array of 

15 v e values (50, 51) includes said array of V r values and an integer number^ 
= V e - V r > 0 of constant values (50)gqual to Os( said number V e of values 
being inversely proportional to the desired pixel size to be obtained for the 

reconstructed image (14). 

2. Method according to claim 1, characterised in that said p 

20 constant values (50) are null valuesJTas = 6| 

3. Method according to claim 1 or 2, characterised in that said p 
values (50) are arranged externally to said array of V r values (51). 

4. Method according to claim 3, characterised in that said p 
values (50) are arranged in a symmetrical way. 

25 5. Method according to claim 3, characterised in that said p 

Q values (50) are arranged in a non-symmetrical way. 

6. Method according to any one of the preceding claims, 
characterised in that the digitized hologram is a rectangular array of V r = 
N r -M r values (51), each value corresponding to a rectangular pixel of sizes 

30 Ax, Ay. 

7. Method according to claim 6, characterised in that the 
hologram reconstructed in the second step is represented by a rectangular 
array of V e = N e M e values, each value corresponding to a rectangular pixel 
of sizes Ac = {XdlNe&x) and Aj/ = (Xd/M E Ay), X being the wavelength of the 

35 wave beam striking the object of which the hologram is recorded, and d 
the distance between the detection device and the object of which the 
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hologram is detected, Ac and being the reconstructed holographic 

image sampling intervals. 

8. Method according to claim 7, characterised in that N e = 

(Ad/ Ax 3 ), Me = (Ad/ Ay 2 ), Ac = Ax, A?? = Ay. 
5 9. Method according to any one of the preceding claims, 

characterised in that, after the second step, if each holographic image 
sampling interval is not equal or less than a certain threshold, the number 
of values /> (50) added to the digitized hologram array is increased and the 
second step is carried out again. 
1(3 10. Method according to claim 9, characterised in that said 

threshold is a function of the signal-to-noise ratio of the holographic image. 

.y.,,^..^'^ 11. Method according to any one of the preceding claims;' 

•p, characterised in that the method is performed for more than one 

holographic images detected at the same time for different wavelength A, 
15 said more than one images being subsequently superposed in order to 
obtain a multi-colour final holographic image (14). 

12. Computer program. characterised in that it comprises code 
means apt to execute, when running on a computer, the method according 
to any one of claims 1 to 1 1 . 
20 13- Memory medium, readable by a computer, storing a 

program, characterised in that the program is the computer program 

according to claim 12. 

14. Apparatus for detection of holographic images, comprising a 
digitized hologram processing unit, characterised in that the processing 
25 unit processes the detected data by using the method according to any 
O one of claims 1 to 1 1 . 
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